The driven Harper model 
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We analyze the driven Harper model, which appears in the problem of tight-binding electrons 
in the Hall configuration (normal to the lattice plane magnetic field plus in- plane electric field). 
The presence of an electric field extends the celebrated Harper model, which is parametrized by 
the Peierls phase, into the driven Harper model, which is additionally parametrized by two Bloch 
frequencies associated with the two components of the electric field. We show that the eigenstates 
of the driven Harper model are either extended or localized, depending on the commensurability of 
Bloch frequencies. This results holds for both rational and irrational values of the Peierls phase. 
In the case of incommensurate Bloch frequencies we provide an estimate for the wave-function 
localization length. 
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1. The Harper Hamiltonian naturally appears in the 
problem of crystal electrons in the presence of a magnetic 
field [l|. It describes the energy spectrum of a tight- 
binding electron in a 2D lattice, whose graphic represen- 
tation is widely known as the Hofstadter butterfly 
The Harper Hamiltonian is a particular case, J x — J y , of 
the more general Aubry- Andre model 



(Hb), = (b l+1 



— J y cos(2iral)bi 



(1) 



where the notations explicitly refer to tight-binding elec- 
trons in a 2D square lattice: J x and J y are the hop- 
ping matrix elements along the primary axes and a is 
the Peierls phase, given by the magnetic flux through 
the elementary cell. If the parameter a in (JXJ) is an irra- 
tional number, the spectrum of H is known to be pure 
point for J y > J x , continuous for J y < J x , and singular 
continuous for J y = J x . Because of this remarkable fea- 
ture the system (QJ also serves as a model of Anderson 
localization in quasicrystals 0, H[ . 

In this letter we discuss the driven Harper model, for 
which the Schrodinger equation for the time-dependent 
quantum amplitudes bi reads 



01 = - — 



e tuJxt bi- 1 )-J y cos(2Tral+u} y t) h 

This model was introduced in our recent publications [5|- 
0] devoted to cold atoms in a 2D optical lattice, subject to 
an artificial magnetic field normal to the lattice plane and 
to an in-plane static (for example, gravitational) force. 
Clearly, this model also describes a tight-binding elec- 
tron in the Hall configuration, so that our results can be 
equally applied to this fundamental solid-state system. 
In this case, the frequencies uj x and lo v are the Bloch fre- 
quencies associated with the components of the electric 
field. 



The energy spectrum and the eigenstates of an elec- 
tron in the Hall configuration crucially depend on the 
commensurability of Bloch frequencies. Namely, for any 
rational ratio u x /u y — r/q (here r, q are co-prime num- 
bers) the energy spectrum is continuous and the eigen- 
states are extended functions [f| Q. It was conjectured 
in Ref. Q that for incommensurate Bloch frequencies the 
energy spectrum is discrete, but the localized eigenfunc- 
tions are characterized by a non-polynomial scaling law, 
with deep implications on the time dynamics of the sys- 
tem. The analysis of the ID system @ presented in 
this paper identifies this scaling law, proves the discrete 
nature of the spectrum for irrational ft — Lo x /uj y and ex- 
plains the interesting dynamical phenomena described in 

2. We begin with a semiclassical analysis of the driven 
Harper model. The classical counterpart of ^ reads 



H c i(t) — —J' x cos(p — uj x t) — J' y cos(x + ojyt) 



(3) 



where p and x are canonically conjugated variables and 
J' = 2iraJ XtV . In fact, it can be shown that Eq. ^ 
follows from (|3j) if p and x are operators obeying the 
commutation relation [x,p] = i2ira, so that the Peierls 
phase a plays the role of an effective Planck constant. 
The system © can be equally studied on the torus (— tt < 
p, x < 7r), on the cylinder (—tt < p < tt, —00 < x < 00), 
and in the plane (—00 < p,x < 00). Considering the last 
case, it is easy to prove that the system is completely 
integrable. Indeed, using the canonical substitution p' = 
p — uj x t and x' — x + uj y t the new Hamiltonian appears 
to be time- independent, 



H', 



-J' x cos(p') — J' cos(x') + uj x x' + uj y p' 



(4) 



and, hence, the right hand side of ((4]) is the global integral 
of the motion. In what follows, however, we shall discuss 
the dynamics of the classical system §5§ on the cylinder, 



2 



that can be compared with the quantum system ([2]) . We 
are interested in the particle motion along the x axis, and 
in particular, in its mean velocity v = lim^oo x(t)/t. 

Although the system © is completely integrable for 
any set of parameters, it has qualitatively different dy- 
namical regimes, depending on whether (3 — u) x /lj v is a 
rational number, and depending on the relative value of 
fi, the system characteristic frequency in the absence of 
driving, 

n = {J' x J' y ) 1/2 ^2^a{J x J y fl\ (5) 

versus to, the geometric sum of the driving frequencies 
to = ^Joj 2 + co y Q. In the high-frequency regime, to O, 

and for irrational (3 any phase trajectory is bounded, 
implying v = 0. The particle can have nonzero mean 
velocity only if /3 is a rational number. This can be 
proved using adiabatic perturbation theory, where one 
distinguishes between the fast variables p',x' and the 
slow variables p, x. We demonstrate this for two par- 
ticular cases: (3 = and (3 = 1. If j3 = the slow 
variable p(t) rs po where po is the initial momentum. 
Then x ~ x + J' x sin(p Q )t, and v = </£ sin(po)- If a 
classical ensemble of particles is uniformly distributed 
over the 'elementary cell' — 7r < p, x < n, we obviously 
obtain ballistic spr eading, whe re the mean-squared dis- 
placement a = \J (x 2 ) — (x) 2 asymptotically follows a 
linear law, <r(t) = At, with A = J' x /^/2. Next, con- 
sider the case (3 = \. As in the former, at zero order 
we have p{t) = p' + io x t = po. However, at first or- 
der, the momentum p(t) = Po + (Jy/^y) cos(x + uj y t) is a 
periodic function of time. Substituting this solution into 
the Hamiltonian equation for the conjugate variable x we 
have x(t) = J' x f Q sm[p(t) — uj x t]6t ~ t, where the propor- 
tionality coefficient can be expressed through the Bcssel 
function Ji(J'/oj v ). This implies that the ensemble of 
particles has a dispersion a(t) = At with A ~ J' x J' y /uj. 
These rates of ballistic spreading are two particular cases 
of a general result, 

A — u;~ (r+9 ~ 1) , w>fi, (6) 

which coincides with the rate of wave-packet spreading 
derived in Ref. @ by using quantum perturbation theory. 

The regime of low-frequency driving, to < Q, is more 
subtle, because here the phase space of ([3]) contains two 
chains of transporting islands, see Fig. [T] Remark that 
these chains exist for both rational and irrational values 
of (3. In a classical ensemble of particles, those with initial 
conditions in the transporting islands move in the nega- 
tive direction at velocity v = uj y , while the others travel 
in the positive direction, and a(t) = At; the values of the 
coefficient A obtained numerically are depicted in Fig. [5] 
for the two values of (3 used in Fig.[TJ As expected, when 
(3 is a rational number, the linear dependence A(w) ~ io, 
valid for to <C fi, leaves place for large ui to the asymp- 
totic dependence ©. For irrational (3 the rate of ballistic 
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FIG. 1: A portion of phase space (stroboscopic map over 
T y = 2n/ujy) of the classical driven Harper model for rational 
(3 = 1/3, left panel, and irrational (3 = (v5 — l)/4 ~ 1/3, 
right panel. The other parameters are J' x y — 2n ■ 0.1545 
and ui = 0.3. Transporting islands are seen as stability 
islands surrounding elliptic points at (x,p) w (0,0) and 
(x,p) w (— 7r, —it). In the case of rational (3 phase trajectories 
are closed on the torus and the stroboscopic map reproduces 
these trajectories. For irrational (3 any trajectory, which does 
not belong a stability island, never repeats itself on the torus 
and appears as a scattered array of points resembling a chaotic 
trajectory. 
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FIG. 2: The rate of ballistic spreading for an ensemble of 
classical particles, versus the driving frequency u>, for rational 
(3 = 1/3 (upper curve, stars) and irrational (3 = (\/5 — l)/4 
(lower curve, open circles). The inset shows the relative size, 
S, of transporting islands (stars) for (3 = (v5 — l)/4 and 
the function (f7| (open circles), where we arbitrary set the 
proportionality coefficient to 4. 

spreading is given by 

A - uS(u) , (7) 

where S(uj) is the relative size of transporting islands 
(1 < S < 0, see inset in Fig. |2J). 

3. We now turn to quantum analysis. An important 
feature of our system is the presence, for w < fi, of chains 
of transporting islands. A system of this kind - the asym- 
metric kicked Harper - was studied in a series of works 
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by Ketzmerick et al. [ltj [TTj], showing that its quan- 
tum dynamics may considerably differ from the classical. 
Here we meet a similar situation, although our system 
has no chaotic component, at difference with that stud- 
ied in [13, [HI]- We have found that, even if the classical 
driven Harper shows a ballistic regime for both ratio- 
nal and irrational /3, the quantum motion of the driven 
Harper is ballistic only for rational /?, while for irrational 
j3 a saturation effect takes place. In other words, for any 
irrational /3, the wave-packet dispersion a(t) follows the 
linear law of the classical system only for finite times, 
being asymptotically bounded from above. The physics 
behind this phenomena is the destructive interference be- 
tween two probability flows going in opposite directions. 
Figure [3] shows the evolution of a localized wave-packet, 
which is initially supported by the central transporting 
island. Tunneling out of this island as well as the op- 
posite process of capture into other islands in the chain 
are clearly seen. The rate of tunneling is defined by the 
ratio between the size of the stability island, S = S(lo), 
and the effective Planck constant fr e ff — 2na. Follow- 
ing analog arguments in Ref. [l(J, we can estimate the 
maximal wave-packet dispersion <j max as follows. Con- 
sider an initially populated transporting island. Due to 
tunneling it is depleted after a time which is exponential 
in S/h e ff 12]. During this time the quantum particle 
is transported at distance ujy/2ira in units of the lattice 
periods. Therefore, 

Gmax ~ — exp ( C^- ] , (8) 

a \ a J 

where C is some constant. It should be mentioned that 
tunneling in and out transporting islands also takes place 
in the case of rational /?. However, for rational /3 the 
interference between two probability flows is constructive 
and a(t) obeys a linear law for all times. 

The observed difference in the wave-packet dynamics 
for rational and irrational j3 indicates a difference in spec- 
tral properties of the evolution operator U. We construct 
this operator by using the substitution bi — » b\ exp(iw x lt), 
which modifies Eq. ([2]) as follows 

ibi = (6{+i + bi-i) - J y cos(27raZ + aj y t)b[ + ujjbt , 

and by integrating Q over the period T y = 2tt /lo v . The 
obtained evolution operator can be approximated to ar- 
bitrary precision by a banded matrix, whose bandwidth 
depends on the parameter J x . It is instructive to con- 
sider the case J x = 0, where the operator U is a diago- 
nal matrix with the elements Uu = exp(— i2ir(3l), which 
is a periodic (aperiodic) function of I for rational (irra- 
tional) (3. [Notice that because of integration over time 
the parameter a does note appear in the last expression. 
This explains or, at least, gives a hint why the driven 
Harper model is insensitive to rationality of the param- 
eter a.] At this point we can draw an analogy with a 
paradigmatic model of quantum chaos - the kicked rotor 
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FIG. 3: Grey tone image of the quantum wave-packet (black 
maximum) versus space I and time t. The left panel is for 
irrational P = (\/5 — l)/4 w 1/3, the right panel has /3 = 1/3. 
The other parameters are J x = Jy = 1, a = 0.1545 and 
lj = 0.45. 
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FIG. 4: The localization length of eigenstates of the evolution 
operator for /3 = (VE - l)/4. The mean participation ratio 
is shown as the function of ui. The inset shows a typical 
eigenstate for ui — 0.6. 

13]. For vanishing kick amplitude, the evolution oper- 
ator of the kicked rotor is also a diagonal matrix with 
matrix elements Uij = exp(— i27r£Z 2 ), which is a periodic 
or aperiodic function of I according to rationality of the 
parameter £. The number theoretic characteristics of £ 
are crucial [14j . since they determine whether the eigen- 
functions are extended or localized. The driven Harper 
model shows the same feature, since the eigenfunctions of 
the evolution operator are localized functions of I if the 
parameter j3 is an irrational number (see inset in Fig. 2]) . 

We studied numerically the localization length as a 



4 



function of the parameter uj. Figure [4] shows the par- 
ticipation ratio P = 1/ J2i \bi\ 4 averaged over 300 eigen- 
states. One clearly sees an exponential increase in the 
localization length with decrease of the driving frequency 
uj, i.e., with increase of the size of transporting islands. 
A similar exponential dep endence is also obtained when 
we vary a at fixed uj [15j. These results confirm the es- 
timate ©, which can be equally used for the maximal 
dispersion (saturation level) and the mean localization 
length. 

4- Conclusions. In the physical problem of tight- 
binding electrons in the Hall configuration, the transport- 
ing states, which are associated with classical islands, are 
responsible for quantum transport of electrons in a direc- 
tion perpendicular to the vector of the electric field, i.e., 
for the Hall current [5Hl6j. It was numerically observed in 
Ref. [6] that for irrational directions of the electric field 
this transport abruptly diminishes when the magnetic 
field (which in the tight-binding approximation defines 
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